Determination of a Wave Function Functional 
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We propose expanding the space of variations in traditional variational calculations for the energy 
by considering the wave function to be a functional of a set of functions X •''/' = rather than 
a function. A constrained search in a subspace over all functions x such that the functional i/)[x] 
satisfies a sum rule or leads to a physical observable is then performed. An upper bound to the 
energy is subsequently obtained by variational minimization. The rigorous construction of such a 
constrained-search-variational wave function functional is demonstrated. 



One of the mostly extensively employed and accurate 
approximation methods in quantum mechanics is the 
variational principle for the energy. Consider a quantum 
mechanical system with Hamiltonian operator H . The 
ground state eigenenergies E and eigenfunctions 5* for 
this system satisfy the Schrodinger equation — E"^ . 
Next define the functional = / 'ip*H'ipdT/ J ^l)*4'dT. 
Searching over all bounded, quadratically integrable 
functions one obtains = E, with ^E" being the 

solution of the Schrodinger equation associated with 
energy E. Since 5I[^] — 0, the functional /['i/'] is sta- 
tionary for ij) = . In practice, an approximate function 
"0 of a particular analytical form is chosen to depend 
upon a number of variable parameters Ci{i — 
A least upper bound to the ground state energy Eq is 
obtained by solving the p equations dl[^j\/ dci = 0, and 
employing the derived set of values of the parameters 
Ci to calculate /[V'] . In application of the variational 
principle, however, the space of variations is limited 
by the choice of form of the function chosen for the 
approximate wave function. For example, if Gaussian or 
Slater-type orbitals are employed, then the variational 
space is limited to such functions. In this paper we pro- 
pose the idea of overcoming this limitation by expanding 
the space over which the variations arc performed. This 
then allows for a greater flexibility for the structure 
of the approximate wave function. We demonstrate 
the idea of expansion of the variational space by example. 



We expand the space of variations by considering the 
approximate wave function to be a functional of the set 
of functions x'- — ^[x]j rather than a function. The 
space of variations is expanded because the functional 
ip[x\ can be adjusted through the function x to repro- 
duce any well behaved function. However, this space 
of variations is still too large for practical purposes, 
and so we consider a subset of this space. In addition 
to the function -0 being of a particular analytical form 
and dependent on the variational parameters Ci , the 
functions x are chosen such that the functional Vix] 
satisfies a constraint. Examples of such constraints 
on the wave function functional f/'fx] are those of 
normalization or the satisfaction of the Fermi-Coulomb 
hole charge sum rule, or the requirement that it lead to 



observables such as the electron density, the diamagnetic 
susceptibility, nuclear magnetic constant or any other 
physical property of interest. A constrained-search over 
all functions x such that V'[x] satisfies a particular con- 
dition is then performed. With the functional '^[x] thus 
determined, the functional J[0[x]] is then minimized 
with respect to the parameters Ci . In this manner 
both a particular system property of interest as well as 
the energy are obtained accurately, the latter being a 
consequence of the variational principle. We refer to this 
way of determining an approximate wave function as the 
constrained-search-variational method. 

As an example of the method we consider its applica- 
tion to the ground state of the Helium atom. In atomic 
units e = h = m = 1, the non-relativistic Hamiltonian is 

H = {-l/2)Vl - {l/2)Vl - Z/n - Z/r2 + l/ri2, (1) 

where ri, r2 are the coordinates of the two electrons, 
ri2 is the distance between them, Z = 2 is the atomic 
number. In terms of the Hylleraas coordinates^ s = 
ri-\-r2, t = ri — r2, u = ri2, we choose the ground state 
wave function functional to be of the general form 



^[x] =$(a,s)[l-/(s,t,u)]. 



(2) 



with $(a,s) 
e-9«(i + gu)[i 



= (a^/7r)ea::p(— as), f{s,t,u) = 
x{q, s, t, u)(l + u/2)], where a and q are 
variational parameters. Any two-electron wave function 
may be expressed in this most general form. 

In our example, we consider x to be a function only 
of the variable s: tp = iplxiljs)]- Then, this wave 
function satisfies the electron-electron cusp condition^, 
which in integral form is ^'(ri, r2, rjv)|ri^r2 = 
^'(r2,r2, ...,rjv)(l -I- ^) + ri2 • C(r2,r3, ...,rjv), where 
ri2 = ri — r2 and C(r2, rs, ta?) is an undetermined 
vector. This is readily seen by Taylor expanding 
'0[x('Z, s)] about u = 0. In a similar manner, the wave 
function ipixilTs)] satisfies the electron-nucleus cusp 
coalescence condition^, which is ^'(r, r2, rAr)|r^o = 
^'(0,r2, ...,rAr)(l - Zru) + r • A(r2, rg, ta,) with 
A(r2, ra, r^r) being an undetermined vector, for 
a = 2. For arbitrary value of the variational pa- 
rameter a, the electron-nucleus coalescence condition 
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is not satisfied. (The differential form of these cusp 
conditions obtained from the integral versions is 



^/dr|r=o — C*lr=Oi where sp.av is the spherical 



average of the wave function about the singularity, 
r = ri — r2, ri and V2 the positions of the two particles, 
and with C = 1/2 for electron-electron coalescence and 
C = — Z for electron- nucleus coalescence.) 

The constraint to be applied to the wave function func- 
tional ■(/'[x] is the normalization condition: 



(3) 



A pictorial representation of the space of variation in 
the above example relative to the standard variational 
method is given in Fig.l. 

The next step is the constrained search over functions 
x(5, s) for which the condition of normalization is satis- 
fied. Substitution of our into Eq.(3) leads to 



f-OQ 

2n^ ds|$(a, 5)1^5(5) = 0, 
Jo 



where 



g{s)^ udu dt{s' -t')[f'{s,u)-2f{s,u)]. (5) 
Jo Ja 

Next consider Eq.(4). If the parameter a were fixed 
at say the value for which the energy is minimized by 
the Hydrogenic pre-factor $(a,s), i.e. a ~ 27/16, then 
there exist many functions g{s) for which the condition 
of Eq.(4) is satisfied. These solutions correspond to the 
subspace C as shown in Fig.l. 

On the other hand, if the parameter a is variable, then 
the only way in which the condition of Eq.(4) can be 
satisfied is if 



g{s) = 0. 



The requirement of Eq. (6) is equivalent to the constrained 
search over the entire subspace C. Substitution of /(s, u) 
into Eq.(6) in turn is equivalent to a quadratic equation 
for the function x{li s)'. 



a-iq, s)xiq, + 2b(q, s)x(g, s) + c{q, s) = 0, (7) 



where 



a{q,s)= / {s^u^ -uy3){l + u/2f{l + qu)^e^^''"du, 
Jo 

(8) 



b{q,s) = -/ {s^u^ -uy3){l + u/2){l + qu) 
Jo 

[e-29"(l + gu)-e-«"]du. 
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FIG. 1: In the figure, the area of box A represents the space 
of all normalized square-integrable functions. Functions such 
as Gaussian or Slater orbitals are represented by lines. The 
space A is also representative of the functionals tl^lxiSjijU)]. 
The subspace area B represents all normalized wave function 
functionals tp[x{s,t)], and its subspace C that of all normal- 
ized functionals i/)[x(s)] • 



(4) c{q,s)^ {s^u^-u^/-i){l+qu)[e-^'i'^{l+qu)-2e-'i'^]du. 
Jo 

(10) 

The integrals for the coefficients a, 6, and c are solvable 
analytically. Solution of the quadratic equation is then 
equivalent to searching over the entire subspace C. 
Subspace C is comprised of only two points. These two 
points corresponding to the two solutions Xiiq^s) and 
X2[q,s) are such that the two wave functions f/'lxi] Etnd 
'0[X2] are normalized. Rigorous upper bounds to the 
ground state energy are then obtained by variational 
minimization of the functional /['(/'[x]] with respect to the 
parameters a and q. The details of the derivation and 
analytical expressions for the wave function functionals 
are to be given elsewhere. 



The ground state energies obtained from the wave 
function functionals V'lxi] and V'[X2] are given in Table 
(6) I together with those due to the energy-minimized 
Hydrogenic pre-factor <I>(s) with a = 27/16, the Hartree- 
Fock(HF) valueQ, and the 'exact' value of Pekeris 4J. 
The functions Xi('?i'5) and X2(9, s) at the respective 
energy minimum are plotted in Fig. 2. We also quote 
energy minimized values determined from the wave func- 
tion functionals ^'[xa] and ^^Ixa] obtained by solution 
of the quadratic equation with a fixed at a = 27/16, 
and q treated as the only variational parameter. (The 
functions (xa, X4) differ minimally from those of (xi, 
X2)-) The satisfaction of the virial theorem and the 
expectation value of the operator W ^ ri + r2 are also 
quoted. 

We note the following points of interest: (i) The 
functions Xi{q-i^) 1 X2{q,s) are very different: Xi{qj^) 
(9) is a positive monotonically decreasing function whereas 
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XiiQjs) is a negative monotonically increasing function. 
They, however, are of about the same magnitude. 
Thus, the functionals i/'ixi] ^^^id '(/'[X2] Siie different 
from each other. Similarly, the functionals V'lxa] ^-nd 
^[X4] differ. (ii) In spite of the wave function functionals 
ipixiilT^)] ^-iid V'[X2('?i s)] being different, the corre- 
sponding energies are essentially equivalent, (iii) Note 
that E[x2] < E[xi] < E[x4] < E[x3]. (iv) The improve- 
ment in the energies noted in (iii) is also reflected in the 
corresponding satisfaction of the virial theorem, (v) As 
the constrained search is over the normalized functionals 
of subspace C, both solutions of the quadratic equation 
Eq.(7) lead to meaningful energies. 

Our results demonstrate the advantage of the concept 
of a wave function functional. (The purpose of the 
paper is not to determine the most accurate wave 
function for the He atom ground state.) The energies 
obtained via the various wave function functionals 
are to our knowledge the most accurate one- and 
two-parameter results in the literature. A 1.93% error 
of the Hydrogenic pre-factor is reduced to errors of 
0.45% — 0.43% for the four wave function functionals. 
The corresponding satisfaction of the virial theorem is 
1.61% — 0.08%. As a further point of comparison, we 
note that our results are superior to those of HF, and 
the results of V'ixi]? 4'[X2] and ipixi] are also superior 
to a 3-parameter wave function calculation(CK) [^. 
( This wave function is of the form of Eq.(2) with 
/(s,u) = e-«"(l + qu)[l - Pe-^f'^l + u/2)] with 
a = 27/16, q = 0.562326, fi = 0.099947, /3 = 0.88066). 
The exact satisfaction of the virial condition by the 
pre-factor is a consequence of scaling, and that of HF 
due to the self-consistent nature of the solution. That 
of Pekeris is a consequence of the accuracy of the wave 
function. The expectation values oi W = ri + r2 are 
of course only expected to be as accurate as that of 
wave function, and not correct to second order as is 
the energy. The values obtained by the functionals 
^[Xi] and V'[X2] are an improvement over those of the 
pre-factor, reducing an error of 4.3% to 2.9%. They 
are also significantly superior to those obtained by the 
3-parameter wave function which is in error by 3.9%. 
Similar improvements are observed for the expectations 
of other single particle operators. On the other hand, 
the high accuracy of the HF value is because in this 
theory single-particle expectation values are correct to 
second order ^ . Our results also demonstrate that by 
expanding the space of variations, more accurate results 
for both the energy as well as other properties can be 
obtained with fewer variational parameters. 

An improvement over the present results can be 
achieved as follows, (i) Expand the space of variations 
by considering % to be a function of the variables s and 
t, or expand the space still further by considering x to be 



FIG. 2: The functions Xi(9iS) and X2(9, s) • 



a function of the variables s, t, and u. Such an expansion 
of the space could require the solution of an integral 
equation for the determination of x- (ii) Replace the 
Hydrogenic pre-factor <I>(q!, s) by the analytical HF wave 
function, (iii) Combine the expansion of the variational 
space with the improvement of the pre-factor. (iv) 
Employ a different analytical form for the wave function 
functional. 

As noted above, it is also possible to search over 
all functions x('?j s) such that the functional ip[x\ 
leads to a physical property of interest. For exam- 
ple, let us consider the expectation of VF = ri + r2- 
< ri + r2 >—< s >= So + As. Here sq is the ex- 
pectation from the pre-factor <i>(s). Assuming As 
known from experiment or some accurate calculation, 
and if a wave function functional V[x] of the form 
in the above calculation is employed, then two dis- 
tinct x's such that < '(/'[x] |s|^[x] >= sq + As can 
be obtained by solution of the quadratic equation 
a{q, s)x{q, sf + 2b{q, s)x((?, s) + [c{q, s) - A] = 0, where 
the constant A — 2As/a*, and where the coefhcients a, b 
and c are the same as in Eq.(8-10). With the functionals 
ip[Xi] and ip[x2] thus determined, the energy could then 
be obtained by minimization of the functional /['(/'[x]] 
with respect to the parameters a and q. In this manner, 
the two wave function functionals would reproduce both 
the size of the atom exactly and the energy accurately. 

For completeness we note that the concept of con- 
strained search in the present work differs from that[3 
within density functional theory (DFT). The idea 
underlying DFT is based on the first Hohenberg-Kohn 
theorem|g according to which the wave function ^ is a 
functional of the ground state density p{r) : "if = '^[p]. 
Thus the energy is a unique functional of p(r) : E — E[p]. 
The in principle constrained search within DFT is as 
follows. One first searches for the infimum of the (H) 
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TABLE I: Rigorous upper bounds to the ground state energy 
of the Hehum atom in (a.u) for different wave functions ( WF) . 
The satisfaction of the virial theorem, and the expectation 
values of W = ri + r2 are also quoted. 
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a = 1.6875 


-2.84766 


2.0000 


1.7778 


i>\x3] 


a = 1.6875 q = 0.581 


-2.89004 


1.9678 


1.7778 




Q = 1.6875 q = 0.180 


-2.89061 


1.9707 


1.7778 




a = 1.6614 q = 0.5333 


-2.89072 


1.9973 


1.8057 




a = 1.6629 g = 0.17049 


-2.89122 


1.9984 


1.8041 


HF 


See Ref. 3 


-2.86168 


2.0000 


1.8545 


CK 


See Ref. 5 


-2.89007 


1.9890 


1.7848 


Pekeris 


See Ref. 4 


-2.90372 


2.0000 


1.8589 



over all antisymmetric , normalized, A^-particle functions 
whose density is p(r). One then searches over all 
such p(r) to obtain the infimum of that expectation. 
The consecutive infima can be shown to be a minimum 
and equal to the ground state energy. In our work, the 
statement that the wave function is a functional of the 
functions x is more general. The functions x ^^re not 
restricted to just being the density or functions of the 
density. And the constrained search over all functions 
X is such that the wave function leads to an arbitrary 
property of interest. The energy is subsequently ob- 
tained by variational minimization. 

In addition to its use in Schrodinger theory, the 
wave function functional ■0[x] ni^-Y ^-Iso be employed 
within Quantal density functional theory (Q-DFT)|s|. 
In Q-DFT, as in traditional Kohn-Sham (KS) DFT|inl|. 
the system of electrons is transformed into one of 
noninteracting fermions such that the same ground state 
density, energy and ionization potential are obtained. 
However, in contrast to KS-DFT which is in terms 
of energy functionals and functional derivatives, the 
Q-DFT framework involves fields whose quantal sources 
are expectations of Hermitian operators taken with 
respect to the wave function. Thus, an approximate 
wave function functional of the form of Eq.(2) can be 
employed in this theory with the Slater determinantal 
pre-factor determined self- consistently. The wave func- 
tion functional ■0[x = f{p)] could also be used. Within 
Q-DFT, the corresponding energy is a rigorous upper 
bound. 

We also mention the work of CoUe and Salvetti(CS) 
[m who suggested a wave function functional of the 
density. The CS wave function is similar to Eq.(2) except 
that the pre-factor is the HF wave function and the 
correlation term is /(ri, =6^*^ [1 — + ^/2)], 

with r = n - ra, R = ri + ra, /3 = q[p"^{Ii)Y'^. 
CS further assumed that the corresponding Dirac 
density matrix was that due to the pre-factor. The 



function x(R) was to be determined by requiring that 
the correction to the HF Dirac density matrix due 
to the correlation factor /(ri,r2) vanishes. Since the 
HF Dirac density matrix cannot be the same as the 
exact interacting system density matrix, this is not 
an exact constraint. Moreover, they did not satisfy 
this condition, and instead approximated the function 
x(R) by x(R-) = + V^/^)- A consequence 

of this was that the resulting wave function was not 
normalizedjlElEi- There is also no discussion in this 
work of the general concept of wave function functionals 
or the idea of constrained search to obtain them. 

In conclusion, we have proposed the idea of expanding 
the space of variations beyond that of standard varia- 
tional calculations by considering the wave function to 
be a functional instead of a function, a functional of 
the functions %. A constrained search is performed over 
the functions x such that the wave function satisfies a 
constraint or leads to a physical observable. A rigorous 
upper bound to the energy is then obtained by varia- 
tional minimization with respect to any parameters in 
the wave function functional. The construction of such a 
constrained-search-variational wave function functional 
for the ground state of the Helium atom where the 
search is over the entire requisite subspace is explicitly 
demonstrated. 
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